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SUMMARY 



Values of the lat eral- st atil i ty derivative C^^, the 

rate of chan£^e of yawing-mcinent coefficient with yav/ing 
angular velocity, contributed 'by the wing, the fuselage, 
and the vertical tail have been det;erinined for a Inid^^^ing 
airiDlane model "037 the f r e c- c 3 ci 11 at i on racthodo 

It was found that the values of C^-^^ contrituted hy 

the vertical tail and by the profile drag of the wing wore 
in good agreement wilh theory. The damping contrihuted 
"by the wing ^';aried as the sauare of the lift coefficient, 
hut the actual values were somewhat lower than those pre- 
dicted "by existing theory.. The value of C-^^ contributed 

by the fuselage appeared to he negligible. 

An empirical formala is presented for obtaining an 
approximate value of C^-^^ for a conventional midwing air- 
plane • 

INTRODUCTION 



*Ii3 calculating the lateral stability of an airplane, 
difficulty is often experienced in estimating values of 
the stability derivative ^n-^> the rate of change of 
yawing-moraent coefficient with yawing angular velocity. Al 
though theoretical methods for obtain ing the value cf ^n-** 
contributed by the vertical tail and the wiag are given in 
references 1, 2, and 3, little recent experimental v/ork 
has been dene to determine values of this derivative. In 
order to provide experimental data on the contributions of 
the wing, the fuselage, and the vertical tail to "^n^^ 
some measurements for a midwing airplane model have been 
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made in the HACA free^flif-ht tunnel. Addit ional ^ measur e- 
monts vjere made for a rectangular y^ring of high^lift sec- 
tion in order to extend the lift coefficients to the high 
values encountered hy full-^scale airplanes. The results 
are presented in the present report* 

A f ree^oscillation method similar to that described ^ 
in reference 4 was used. The values of Cnp were direct- 
ly determined from the dampiag of free^-yawing oscilla- 
tions, which were obtained with the models mounted on a 
strut that permitted freedom only in yaw. 



CO 
-0 , 



SYMBOLS 



AO 
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rate of change of yawi ng-moment coefficient with 
yawing angular velocity per unit of rb/sV 

rate of change of yawing--moment coefficient v/ith 
angle of sideslip (dCj^/36) 

lift coeiiTcient (L/qS) 

lift coefficient of wing alone 

increment of lift coefficient due to flap 

prcfile-^drag coefficient (D^/qS) 

profile-drag coefficient of wing alone 

increment of profile-drag coefficient due to flap 

yav/ing-moment coefficient (lT/o;^bS) 

yawing mom.entj foot-pounds 

rate of change of aerodynamic yawing moment with 
yav/ing angular velocity (oIT/dr) 

rato of change of frictional yawing moment with yaw- 
ing angular velocity [(di^^/dr)^] 

rate of change of aerodynamic yawing moment with an- 
gle of yav/ (dN/dv) 



rate of change of roBtoring moraent of torsion spring 
wiGh angle oi yaw 

lift, J pounds 

profile drags pounds 

dynamic; pressure, pound3 per square foot (^gp^ j 
wing area, s^.^uare feet 

yawing angular velocity^ radians per second 

wing spanj> feet 

flp.p span, feet 

airspeed, foot per second 

air density, slugs per cul^ic focb 

angle of ya.Wj radia.ns 

maximum amplitude of yawing oscillation at zero 
time, ra.dians 

maximum amplitude of yawing oscillation at time t, 
radians 

angle of sideslip, radians C-\|/) 

total logarithmic decrement or damping factor 

logarithmic decrement due to friction 

time; seconds 

period of yawing oscillation, seconds 
aspect ratio 

taper ratio (ratio of tip chord to root chord) 

distance from center of gravity to rudder hinge line, 
feet 

yawing moment of inertia, slug-feet square 
, k-j . kf constants 
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METHOD 



The equation of molicn of a system heving freedom only 
in yaw can be expressed, to a close first approximation, as 

I, ^'4 - fM. iT,.^ ^ - (n^ . ^^ = 0 • (1^ 

^ db'" \ ^ dt V V y 

The yawing motion of the system represented by this 
equation can he expref^sed by an equation of the form 

\1/ = e " ^ (A sin b t + B cos d t ) 

v/hich represents a damped harnonic oscillation v/here the 
ratio nf' the maximum amplioudes of euccessive oscillations 
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max4. -ai; 
e 



^max^ 



The value of a, the losarithriic decrement or th3 damping 
factor, can be determined frcn. the experimentally recorded 
angles of ya^;^ against bime by. means of this relationship 
which t \>rhen tra^nspcsodj giv6r3 

a ^- V ^ /■ 



The damping derivative expressed in terms of the damp- 
ing factor is 



f 

and the damping derivative due to friction is 



Combining equ.ations (5) and (4y gives 

ITp = -2l;7; (a - af) 
or,, in ncndimonsiohal form. 



C._ = . 1I4 (a -. a,) ' 



' n Q -K 
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The period of the vav/ing oscillation expressed in 
terms of the coefficients of equation (i) is 



S T = ^ = (6) 



The effect 01 friction on the period is negligitlo. At 
zero airspeed, when Nj. and IT^, "become zero^ equation 
(6) reduces to 



or 



Z 
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By substituting i a equation (?) the \'-alue of T" at 
zero airspeed, the ya'vving moment of inortis. '^^ 
obtained for use in equation (o). 

It should he noted that the restoring moment of the 
torsion sprinf^ k affects the period of the oscillation 
(equation (6)) hut does not affect the damping (equation 
(3)), It is possible^ therefore, to adjust the period to 
any desired value nifchout affecting the measurement of 

^nr • 

APPARATUS AFD PROCSDUSiiJ 



The investigation v/as carried out in the NACA free- 
flight tunnel with the apparatus shown in figure 1, The 
upper portion of the struo to which the model is attached 
is mounted m hall hearings and is free to rotate within 
the fixed hase. The model is therefore free to yaw hut is 
restrained in roll and pibch. The movable portion of the 
strut is hinged to permit adjustments in the angle of at- 
tack of the model being tested. 

A torsion spring connecting the fixed and movable 
portions of the strut provides the additional restoring mo- 
ment necessary for obtaining short-period yawing oscilla- 
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tions. Tl is iirportant that the period of t he ' o scilla- 
tions "bo fairly slooi't to injure a ^-/el 1 -de i i ne d oscillation 
envelope and therefore to permit an accurate measurement 
of dampings 

The airplane model used in tho" tests is shown in fig- 
ures 1 and 2. Th3 \fing of the niodol had an aspect ratio 
of 6„7 and a taper ratio of 0.40^ and was equipped v/ith 
partial-span split: flaps deflected 60^ • Two vertical 
tails, shown in fi^-ure 2, were used on the model ♦ The 
model v;as mounted on the strut with its center of gravity 
on the axis of rotation. 

''fho rectangular wing used in the investigation had an 
aspect ratio of 6 and an ITACA 103 airfoil section. This 
airfoil section wa'i used tecauce of its high maximum.^ 1 i ft 
coefficient at the low Ileynolds numbers of the f r eo-- f 1 i ght - 
tunnel tests. 5'or Bome of the. tests the wing was fitted 
with a split flap ■ 20 percent of tho wing chord and 60 per- 
cent of the wing span. 

The airplane model was tested at-dynamic pressures of 
lc9 and 4. r pounds per square foot* l^o aiopr e ciahl e change 
^'np was noted v/ith variation in dynaraic pressure. 

The rectangular wing was tested only at a dynamic pressure 
of lo9 pounds per souare foot "because of excessive vibra- 
tion of the wing at higher values of dynamic pressure. 

The testing procedure conr.istid simply in yawing the 
model approximately 10^, releasing it, and recording the 
resulting, oscillations v/ith a mo i on-p i ctur e camera mounts 
ed o.n top of che tunnel. 

The , friction of the system was determined ' from tests 
at zero airspeed with t?ie models replaced "by flat load 
weighus on long rods. These weights were adjusted to sim- 
ulate the yai/ing moments of inertia of the models and were 
alined x-'ith the plane of rotation to give negligible air 
damping^ In tests of the airplane model at 2:oro airspeed 
with vertical tail removed, essentially the same damping 
was ohtained as in the friction test^^o It appeared, there-- 
fore, that a tail-off run at' zero airspood could be satis- 
factorily used to replace the cpecial friction tests with 
lead weights. 

The peaks of the oscillations recorded hy the camera 
were read from the film record and plotted against time. 
The natural logarithms of the faired peaks were then plot- 
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ted against time and the slopes of the resulting straight 
linos were graphical x epro sen t a t i on n of the 1 og^-^i thmi o 
decrements a and a^-.. 7he numerical values for a and 

a-f were determined irom the slopes hy eqi:.ali on (2) and 

CO these values were sul:) s 1 1 tut ed in equation (5) to obtain 

I ^ n V, • 

P. ^ 

Lift and drag coefficients and yawing-inoment coeffi- 
cients due to sideslip were determined "by tests on the 
si x-.-'Con3ponenk halance in the ttmnel for use in correlat- 
ing the measured values of Oji. , with the nheoretical de- 
rivative Sc 

TflSOR^lTTCAL DAMPING DJIHIVAT ITS S 



The value of Cjj^, for a complete airplane may he as- 
sumed bo he made up of directly additive c jnfcrihut ions of 
the vertical tail, vringj and fuselage, if interference ef- 
fects arc neglected; that is^ 

Cn^ - ^-^-(tali: '^^^^r(wing) ^'^^r( fuselage ) 

It can "be shown ohat the cont r ihut i or. of the vertical 
tail 1 3 

^^n / 2 ~ AC^. . . (S) 

^r(tail) h -^'3 (tail) 

"Fov a wing without flaps 

^^r ( ^^rm g) o a. Q . J. L, 
Simple integration for IIq yields 

Values for are given in figure 13 of reference 2s 

vhicli may he represented "by the Gquaticri 

K, = ^0.031 (1 ' (10) 

V 13 2-5 /' 

The value -C.033 is for a rectangular v.-ing of aspect 
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ratio 6.0. G-lauert, in reference 1, gives a value of 
-0.024 for this condition. 

For a win^: v/ibh partial-span flaps extended, the 
profile-drag term K^Cj)^ "beccmes 



Kq^D - -o^B^ - K^l^Cj, (11) 

^ ^f . 



where 



'A (12) 



and the induced-drag term KiCl^ takes the form 



(13) 



Values for and K3 are given in figures 12 and 

13 of reference 2, out the value for is not given in 

this reference and Is apparently not c^vailahle from other 
sources. Inasmuch as AC^ , . k for the flaps-extended 

condition cannot ho computed v/ithout the value of , it 

aTjpears desirable to prepare additional charts for this 
factor. 

Calculations of reference 5 indicate that the theoret- 
ical value of AC^ . .i.,-.p^ ^ero for fuselages that 

jj \ .L u s e -L r.^ e / 

are ellipsoidal in shape. 

HB3ULTS AjD DISOU'SSIOIT 
Contri Dution of Vortical Tail to 

Valuer of ^r... the complete model ivith partial- 

span fia-ps extended are given in figure 3 as a function of 
vertical --tail siz.e. Values of ^Cnr^(t,ail) ^^'^ obtained 
directly from the data in figure 3 by subtracting the 
value of Cr with tail off from the values of v/ith 

tail on. The line drawn on the figure was computed from 

eouation (8) and was based on bhe measured values of 

AC" / N given in the table on figure 3. The agreement 

k t ail ; 
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"bet^'/een the test points and this computed line is an indi- 
cation that the cont rihut i on of the vertical tail to C^^^ 

for a. midx'/ing a-irplane can "be computed vith reasonahle ac- 
curacy from the theoretical relation given in eauation 
CO (s) . For high-- or low-wmg airplanes a correction factor 

might "bo necessary for this relation "because the sidewash 
at the tail, which varies with i\^ing position, causes dif- 
ferent changes in AC^^ , . and ACn • v. 

g(tail) r(tail) 

Contrioution of ITuselage to C 

n p 

The data of figure 4 show the variation of C-^^ with 

lift coefficient for the fuselage-^wing conhination and for 
the rectangular wing with partial-span flaps extended. 
The value of C^^^ for the fuselage and v;ing varied from 

-0.014 to ^-0.028 over the lift ra.nge covered in the air- 
plane model tests. 

A comparison of the 0^.^., values for the fuselage- 
wing com-hinaticn v;ith the values for the rectangular wing 
v:ith flarjs extended (fig. 4) indicates that the fuselage 
had a negligihle effect on ^'n^ • Although it appears from 

a direct comparison cf the data that the fuselage slight- 
ly reduced ^n^f this apparent reduction v/as probably 

caused "by the difference in plan form and "by the greater 
profile drag of the rectangular wing. Other recent tests, 
the results of vrhich ar e 'iinpuhl i she d , have indicated val- 
ues of AC^.^ , , ranging from -0,003 to --O^OOS. It 
rv fuselage/ 

appears, therefore, that the fuselage contribution to 

Hp 

is normally small enough to he neglected. 

Oontrihut:on of V/ing to 

Variation of ACv. , , v with lift coefficient.- The 
^ r V wing; 

data cf figure 5 shov/ that ACv. . . ^ for the wing with 

•^^r(wing> 

flaps retracted varied as the sq^uare of the lift coeffi- 
cient, as predicted hy theory, hut that the value of K.^^ 
was smaller than the value predicted hy either reference 1 
or reference 2. The experimentally determined value of 

for the rectangular wing was -0.020; whereas reference 
1 predicted a value of -0.C24 and reference 2, a value of 
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-0.031. It appears that the value of -0.031 given liy ref- 
erence 2 and used in equation (lO) is too large and should 
be replaced hy -0.020* 

Tno variation of 0^,^ wSth lift coefficient for the 

v/ing with par bial-^span flaps er.tended (fi^.. 4) differed 
from the variation with flaps retracted in that the mini- 
mura value of C^, v:as ohtaincd at a small positive lift 

coefficient rather than at zero lift.. This rosulT:. v;hich 
is also indicated by equation (l3), is due to the fact 
that at zero lift the center flapped section is developing 
positive lift, the tip section is developing negative lift, 
and hoth are cent r1. hut ing to <^nj, > Inasmuch as no calcU'- 
latcd value for the constant B'^ vras available, no corre- 
lation of the theoretical and exper ircental variation of 

AG,. , V with lilt coefficient could he made for the 

^^r ( wing ) 

f la.ps-ext ended condit i on . 

Variation of . , ^ with profile drag,- The 
^ r \ \) 1 n g J 

value of for the v/^nir ./ith flaps retracted at zero 

lift was about -0.007, as shovrn in figure 5. The profile- 
drag coefficient Cy^^ for the win-^, as measured on the 

balance in tne tunnel for the same condition, was 0.084. 
J'rora these two values, h^ is found to be 

-0,29. Equation (9) yields 0,33 as the theoretical value 
for for a rectangular wing. It aopears "chat the cal- 

culated and the experimentally determined values of 
are in fairly good agreement. 

With the partial-span flaps deflected on the rectangu- 
lar wing, the value of Cp_ due to profile drag can be ob- 
tained from the value of Cv. at the lift coo:^''f icient 

*r 

given by the flap. Tor the wing tested, the flap gave an 
increment of lift coefficient of 0.60, ?rcm figure 4 at a 
lift coefficient of 0.60 the value of 0^.^ was -0.017. 
Combining equations (ll) and (13) and eliminating terms 

containing Ov , because Cr = 0 at Oj^ = 0.60, gives 
•^^■^f -^w 

■ r 0^^ o J - 



11 



The value ^-rS'T) was --0.007 from the v;ing-alcne test; 
^Ct) was 0.080 iroP force tests: K.-. v/as -<0.072 from 

O 'I 



equation (13); and vas -0,0092 from reference 2. Then, 

for a value of ACl^. of 0,60» 

r 2I ■ 

Oj. - -0.007 -h (^0.072X0.080) + j -0 . 009 3x ( 0 . 50 ) I 

L . J • 



= -0,007 -0.006 -0.00 
.-.0.015 

This result is in frood a^^reement with the measured value of 
-0a017. The nai^ni tude of all of these factors is sinall, 
however I, compared vith the contribution of the tail surface. 

Determination of Cj.^ for Complete Model 

The folloTving empirical formula, which was derived 
from test results, should eive a fair approximation of the 
value of for a conventional midwing airplane; 



coNCLUDiira remapks 

The free-oscillation method of determining damping in 
yaw is considered vcTy satisfactory in that it i::rovided 
reasonahly accurate results quickly and easily. The follow- 
ing conclusions were drajn from the results of the free- 
oscillabion tests of the midv/ing airplane model and the rec- 
tangular V73ng model: 



1* The experimental values fov- the yawing moment due 
to yavjing c on b r i "bu 0 e d "by the vertical tall v/ere in good 
agreement v;ith the calculated values. 

2c The values of the yawing momenb duo to yawing 
contrihuted by the wing varied as the square of the lift 
coefficient "but were lower than those predicted hy theory. 
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3. The value of the yawing moment due to "yawing con- 
tributed "by the profile drag of the wirxg was approximately 
the same as the thooretical value* 

4^ The coritr ihution of the fusela^^e to the yawing 
moment due to yawing was negligihlo compared with the 
value for the complete model, 

5o The tost resiiults indicated that a fair approxima- 
tion ^f the value of the yawing moraont due to yawing for a 
conventional midwing airplane could he obtained from an 
empirical formula „ 

Langley Memorial Aeronautical I-ahorat cry ^ 

National Advisory Gommitt;ee for AGronaatic3» 
Lan^-^ley I'ield, Vac 



R3F£1R3¥0BS 



1. Glauert., H.s Calculation of the Rotary Derivatives 

Due to Yawing for a I'onoplane Wing. S. & NOc 
866, British AcR.G,, 1933. 

2. Pearson, Henry A., and Jones, Robert T.; Theoretical 

Stability and Control Characteristics of Wings with 
Various Amoun-cs of Taper and Twist. Hep. No. 635, 
ITACA, 1938, 

3. Zimmerman, Charles H.? An Analysis of Lateral Stabil- 

ity in Power-Oif 5'iight with Charts for Use in De- 
signo :Rep, ITo • 599, lUCA , 193?. 

4. Dramwell, 3^. H., and Relf, '? . i Experiments on Model 

of Complete Aeroplanes, Part IV - Determination of 
the Pitching Moment Due to Pitching for a Model^ Bi^ 
plane at Various Inclinations to the Wind* R. & M. 
No, 111, British A. C .A . , 1914. 

5. Kunk, Max II. i Fundamentals of J'luid Dynamics for Air- 

craft Designers^. The Ronald Press Cov, 1929, 



L~3C7 




\ 



"...■Partial- -pan split flap 
O.oO sps-n, 0.25 chord 



Fi^;iire 2,- Silo elevation of rnodol rsr.d in lamping tests in IvACA f reo-f li^-ht tannel. 
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Fig. 3 




_ Variation of lanipi-ag in jav/ vn th vertical- tail efiQctive- 
ness, Midwing airplane nxodelj jflaps extended^ = 1,0. 
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Figure 4.- Variation of iainping in yav/ with lift coefficient. 
iTlapo exteniocl; vertical tail reinovel. 
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Sz-rperi-nental values 
-(0. 0074-0, 0?0 Cl^) ' 
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^(0.0074-0.031 CL^)(ref. P) 
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Variation of da^r^ping in yaw \;itli lift coefficient. 
Hectan^^lar v/ingj flaps retracted. 



